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I.  Introduction 

A  variety  of  Air  Force  applications  require  highly  detailed  information  about  atoms, 
molecules,  and  their  interactions.  This  information  is  necessary  in  problems  ranging  from  chem¬ 
ical  laser  development,  to  the  detection  and  identification  of  rocket  plumes,  to  non-linear  optics, 
electron  beam  technology,  and  high  density  and  energy  fuels. 

The  crucial  component  needed  to  understand  molecular  reactions  is  the  potential  energy 
surfaces  (PES)  that  serve  to  describe  the  attractions  among  the  atoms  and  molecules.  However, 
such  information  is  not  easy  to  obtain.  In  many  cases,  the  most  direct  approach  to  obtaining 
accurate  potential  surfaces  for  molecules,  and  detailed  information  about  their  excited  states, 
vibrational  spectra,  and  a  wealth  of  other  quantities,  is  high  level  ab  initio  solultions  of  the 
Schrddinger  equation. 

However,  more  so  than  in  most  other  areas,  the  ability  to  provide  reliable  quantum  mechani¬ 
cal  results  for  increasingly  large  molecules  depends  critically  on  improved  method  development. 
Whereas  supercomputers  can  enable  us  to  make  much  larger  computations  with  old  methods, 
the  simultaneous  development  of  new  methods  can  increase  computational  capability  by  further 
orders  of  magnitude.  In  this  regard,  many-body  perturbation  theory  (MBPT)  and  its  infinite-order 
extensions  termed  coupled-cluster  (CC)  methods  offer  a  number  of  attractive  features  that  the 
more  traditional  configuration  interaction  approaches  lack. 

Under  AFOSR  support,  we  have  established  these  CCyMBPT  theories  as  being  among  the 
most  accurate  available,  and  have  developed  very  efficient  and  generally  applicable  computer 
programs  to  perform  CC/MBPT  calculations.  Also,  we  have  employed  these  methods  for  the 
first  time  in  large-scale  ab  initio  calculations  of  potential  energy  surfaces.  The  successes  of  our 
original  work  in  this  effort  have  been  substantial  (see  previous  AFOSR  reports). 

Following  a  statement  of  research  objectives,  we  review  the  recent  scientific  advances  we 
have  made  under  AFOSR  support 

II.  Research  Objectives 

The  overall  objective  of  our  continuing  research  program  for  AFOSR  has  included  the 
following: 

1.  The  development  of  new,  more  accurate  and  more  efficient  ab  initio  quantum  mechanical 
methods  based  upon  coupled-cluster  (CC)  and  many-body  penurbation  theory  (MBPT)  for 
determining  molecular  properties  and  potential  energy  surfaces  for  molecular  interactions. 

2.  The  implementation  of  these  methods  into  highly  efficient,  transportable  computer  programs 
(ACES  II),  to  enable  computations  on  molecules  to  be  made  on  an  almost  routine  basis,  for 
a  wide  variety  of  different  properties. 

3.  The  application  of  these  techniques  to  a  variety  of  problems  that  are  of  interest  to  AFOSR. 
and  that  serve  to  establish  the  range  of  accuracy  for  CC/MBPT  methods. 

The  underlying  goal  of  our  effort  is  highly  accurate  first  principle  quantum  mechanical 
predictions  of  molecular  properties.  Such  tools  are  already  making  possible  the  determination 
of  dissociation  energies,  moments,  polarizabilities,  excitation  energies,  and  vibrational  and 
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rotational  spectra  for  molecules  to  an  accuracy  that  is  predictive  and  comparable  to  experiment. 
Furthermore,  for  transient  molecules  as  occur  in  rocket  plumes,  combustion,  flame  chemistry,  and 
interstellar  space  that  are  generally  difficult  to  study  experimentally,  theoretical  calculations  will 
provide  otherwise  inaccessible  information.  In  addition,  theoretical  work  provides  the  underlying 
framework  to  understand  innumerable  important  molecular  phenomena  and  to  suggest  previously 
unanticipated  solutions  to  defense  related  problems.  All  of  the  new  methods  which  we  develop 
are  being  incorporated  into  our  Advanced  Concepts  in  Elecu-onic  Structure  11  (ACES  II)  program 
system.  The  program  has  been  established  at  Air  Force  installations  for  use  by  other  Air  Force 
researchers. 

III.  Notable  Accomplishments  This  Past  Year 

A.  A  primary  achievement  of  this  past  year  is  the  implementation  and  application  of  the  full 
CCSDTQ  method  for  several  molecules."**  Such  a  method,  which  includes  the  connected 
T4  cluster,  in  addition  to  Ti,  T2  and  T3  is  exact  for  four-electron  systems,  and  through 
four-electron  clusters  in  many  electron  molecules;  and  is  correct  through  sixth-order  in 
perturbation  theory.  We  introduced  fully  factorized  equations  using  the  recursive  generation 
of  intermediate  approach  we  recently  developed.  Hence,  despite  the  highly  non-linear  nature 
of  the  CCSD(T)  equations,  all  operations  only  involve  multiplication  of  an  integral  with  a 
single  t  amplitude.  Compared  to  full  Cl,  the  average  errors  for  CCSDTQ  are  14  /zhartrees 
at  equilibriuum  geometries  and  26  ^ihartrees  at  twice  R*. 

B.  Another  notable  achievement  this  year  was  our  report  of  the  first,  open-shell  MBPT(2) 
analytical  second  derivative  (Hessian)  method.^  For  unrestricted  Hartree-Fock  (UHF)  ref¬ 
erences,  we  developed  a  general  purpose  program  and  added  it  to  our  ACES  II  program 
system.  We  report  results  for  force  constants  and  vibrational  frequencies  for  the  interesting 
CF2  molecule  in  its  triplet  state. 

C.  The  ketenyl  radical  HC=C=0  plays  a  role  in  a  number  of  important  reactions,  such  as  those 
involving  acetylene  and  oxygen  and  it  occurs  in  interstellar  space.  It  exhibits  a  Renner-Teller 
splitting,  causing  it  to  be  bent  in  its  ground  state,  but  its  detailed  structure  and  vibrational 
frequencies  are  unknown.  We  provide  predictions^  using  our  recently  developed  analytical 
gradient  CC  methods.  As  the  problem  is  complicated  by  symmetry  breaking,  careful  use  of 
symmetry  adapted  coordinates  are  required  to  obtain  second  derivatives. 

D.  Another  advance  this  year  was  the  first  development  of  open-shell  analytical  coupled  cluster 
gradients  for  methods  that  include  triple  excitations.®  We  developed  the  equations  and 
implemented  into  ACES  II  analytical  derivative  methods  for  MBPT(4),  CCSD+T(CCSD), 
CCSD(T)  and  (JCISDCT)  with  a  UHF  open-shell  reference.  The  primary  application  was 
to  the  important  atmospheric  radical,  HOO.  Structures  and  vibrational  frequencies  were 
obtained  at  several  levels  of  approximation. 

E.  In  a  similar  study,  we  generalized  the  analytical  gradient  theory  for  triple  excitation  methods 
to  use  a  restricted  open-shell  Hartree-Fock  (ROHF)  reference  function.^  Several  examples 
including  comparison  with  full  Cl  are  reported. 

References  in  this  section  refer  to  the  AFOSR  publication  list  on  p.  5. 
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F.  In  our  continuing  efforts  toward  defining  the  best  non-iterative  higher-order  correla¬ 
tion  corrections  we  studied  several  examples  where  full  Cl  results  are  available  (BH, 
HF,  H2O).  We  used  the  non-iterative  triples  model  CCSD+T(CCSD),  triples  and  quads, 
CCSD+TQ*(CCSD),  and  the  CCSDT+Q(CCSDT),  non-iterative  quads  augmented  CCSDT 
method,  plus  several  new  methods.’®  Even  in  regions  of  approaching  severe  degeneracy, 
where  a  multi-reference  method  would  be  strongly  preferred,  the  higher-level  single  ref¬ 
erence,  non-iterative  methods  CCSD+TQ*(CCSD)  and  CCSDT+Q(CCSDT)  were  found  to 
perform  exceptionally  well,  comparing  favorably  with  reference  full  Cl  results. 

G.  Like  the  UHF-MBPT(2)  open-shell  second  derivative  method,  another  substantial  achieve¬ 
ment  was  deriving  and  implementing  analytical  MBPT(2)  second  derivatives  for  the  restricted 
open-shell  Hartree  Fock  (ROHF)  based  MBPT  method.^  Unlike  UHF,  ROHF  provides  a 
wavefunction  which  is  a  spin  eigenfunction.  This  makes  it  possible  to  get  much  improved 
results  for  cases  where  UHF  has  a  large  degree  of  spin  contamination.  Our  illustrative  ex¬ 
ample  was  for  the  methylimidogen  (CH2N)  radical,  which  is  strongly  spin  contaminated. 
ROHF-MBPT(2)  gives  much  better  vibrational  frequencies  compared  to  UHF-MBPT(2). 
As  ROHF  MBPT  requires  the  development  of  a  new  MBPT  method,  analytical  second 
derivatives  represent  a  major  accomplishment  and  its  presence  in  ACES  II  guarantees  wide 
application. 

H.  Finally,  we  address  some  formal  questions  concerning  the  nature  of  energy  functionals  and 
the  appropriateness  of  various  alternative  ansiiize  in  CC  theory.^  ”  We  consider  unitary 
(UCC),  expectation  value  (XCC),  and  the  extended  (ECC)  method,  and  their  satisfaction  of 
several  desiderata.  Only  ECC  satisfies  all  requirements,  but  it  is  prohibitively  difficult  to 
apply.  We  find  that  standard  CC  methods  are  difficult  to  top  for  realistic  applications. 

All  of  the  above  are  important,  timely  developments  in  the  field  that  provide  solutions  to 

molecular  problems  that  in  many  cases,  could  not  be  previously  addressed. 
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